Exciton Transfer Integrals Between Polymer Chains 



William Barford* 
Physical and Theoretical Chemistry Laboratory, University of Oxford. 
South Parks Road, Oxford, 0X1 3QZ, United Kingdom 

The line-dipole approximation for the evaluation of the exciton transfer integral, J, between 
J"""""- , conjugated polymer chains is rigorously justified. Using this approximation, as well as the plane- 

wave approximation for the exciton center-of-mass wavefunction, it is shown analytically that J ~ L 
when the chain lengths are smaller than the separation between them, or J ~ L _1 when the chain 
lengths are larger than their separation, where L is the polymer length. Scaling relations are also 
5^ . obtained numerically for the more realistic standing-wave approximation for the exciton center-of- 

mass wavefunction, where it is found that for chain lengths larger than their separation J ~ L 
\ or J ~ I/ -2 , for parallel or collinear chains, respectively. These results have important implications 

for the photo-physics of conjugated polymers and self-assembled molecular systems, as the Davydov 
\Q ' splitting in aggregates and the Forster transfer rate for exciton migration decreases with chain 

lengths larger than their separation. This latter result has obvious deleterious consequences for the 
performance of polymer photovoltaic devices. 

I . I. INTRODUCTION 

The transport of excitons - and hence energy - through molecular materials is determined by the exciton transfer 



<Z3 



integral between neighboring conjugated chromophores. Typically, the length of these chromophores are larger than 
the distances between them, and thus the usual expression for the transfer integrals, derived by invoking the dipolc 
approximation, is invalid. Indeed, the use of the dipole approximation implies that the exciton transfer integral 
between equivalent neighboring chromophores scales linearly as the chromophore length. This prediction contradicts 
the computations of the transfer integrals using various computational techniques, which predict that the transfer 
integrals vanish in the asymptotic limitji^^ - a result in agreement with the decreasing Davydov splitting as a 
O : function of chromophore length observed in self-assembled molecular aggregates^. 

The exciton transfer integral, J mn , between two polymer chains (or linear conjugated chromophores) m and n is 
defined by£ 

1> - 

(TV Jmn= Vij m (GS\Ni\EX) m n (EX|JV J |GS) n , (1) 

m : 

(N ' 

where |GS) m and |EX) m are the ground and excited states of polymer m. is the Coulomb interaction between 
electrons in orbitals cj>i and <fij, and N{ is the number operator for electrons in orbital fy. m (GS| -ZVi|EX) m is thus 
the transition density for orbital i on polymer m. When |r< — r ^ | is large compared to the interatomic spacing the 
Coulomb potential is 

J: v ->'wk\- (2) 

O ' where rj is the coordinate of atom i. 

In practice, the exact computation of exciton transfer integrals using transition densities is computationally expen- 
. £h ! sive. A convenient approximation is the so-called line-dipole approximation 8 '^, which provides a physically intuitive, 
yet accurate description for exciton transfer when the chromophore separation is large enough (typically, three or four 
$H ' times the monomer length). 

In this paper we use the exciton model to explicitly derive the line-dipole approximation. Next, using the line- 
dipole approximation, analytical expressions for the exciton transfer integral between parallel and collinear conjugated 
polymers are derived, and the scaling with polymer length is determined. We conclude by discussing the implications 
of these results for exciton transfer via the Forster mechanism. 

Before deriving the line-dipole approximation, however, we review the point-dipole approximation for the evaluation 
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of J mn- We define r, and ?j as the site coordinates relative to the center-of-mass of their respective molecules, 

Y'L I"i Rm; 

f j = V J - R n> ( 3 ) 

where R m and R„ are the center-of-mass coordinates of molecules m and n, respectively. Then, if 

|fj - rj\ « |R m - R„| = |R„ m | (4) 
we may perform the dipole approximation and write, 



| r . _ r .1 ~ If? 13 IT? 15 ' ^ 



L ^ J. i 

Finally, substituting into Eq. ((!} gives, 

Jrnn — ^rnn J win ' (^) 

where 

o _ [ m (GS|/i m |EX) m ] [ n (EX|ft n |GS) n ] 

mn i-R 13 ' *> ' 

l^mn | 

and 

ft>mn Y m • v n 3(R mn • r m )(R mn • r n ), (8) 

is an orientational factor. r m and R„ m are the unit vector parallels to p. m and R m „, respectively j± m is the electronic 
dipole operator for molecule m, defined by 

Am = e ^ f jNj, (9) 

and 

m (GS|£ m |EX) m = Mr (10) 

is the transition dipole moment of polymer m with N repeat units. 

In practice the condition Eq. Q is far too severe for conjugated polymers in the solid state, and Eq. (J7J) is not 
applicable. However, the much weaker condition that, 

d«|R ron |, (11) 

where d is the monomer size is often satisfied, and under this condition the line-dipole approximation becomes valid. 



II. THE LINE DIPOLE APPROXIMATION 



In this section the line-dipole approximation will be justified for a simplified model of polymers, namely a chain of 
dimers (or double bonds) connected by single bonds, as illustrated in Fig. [TJ This model is of course applicable to 
polyacetylene. However, more generally it is also applicable if the dimer represents a monomer. 

To derive the line-dipole approximation it is also necessary to introduce a model for excitons in conjugated polymers. 
In the weak-coupling limit, defined by the electronic band width being greater than the Coulomb interaction, the 
excited states of conjugated polymers are Mott-Wannier excitons^iS, defined by 

\EX) = Y,1>(r)*(R)\R,r), (12) 

r,R 

where \R,r) is an electron-hole basis state, 

l^' 7 *) = ~?K ( C fl+r/2,T C fl.-r/2,T + C R+r /2,l C R-r /2,l) l GS )- ( 13 ) 
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c R+r/2 ct creat es an electron with spin a in the conduction band Wannier orbital at R + r/2, while c v R _ r j 2 a destroys 
an electron with spin a in the valence band Wannier orbital at R — r/2, and |GS) is the ground state. ip(r) is the 
'hydrogenic' wavefunction for the particle- hole pair, where r is the relative coordinate. ^(R) is the center-of-mass 
envelope wavefunction, where R is the center-of-mass coordinate. 

To a good approximation it can be shown 7 - that the transition densities satisfy, 

(GS|^|EX) = V(R U ) (^Pj (14) 

We note that the transition densities are modulated by the center-of-mass wavefunction and alternate in sign, as 
illustrated for the lowest excited exciton in Fig. [T] 
Using Eq. (|14[) . Eq. fl} now becomes, 

J mn = ^^^•*( J R, i )*(^.)(-l) i + 1 (-l)i+ 1 , (15) 



To simplify this expression it is convenient to partition the sum over sites, i, as a sum over unit cells, t, and a sum 
over sites within a unit cell: 

Jmn = (yf) £ £ J E E ■ (16) 

Then, if the unit cell size, d, satisfies, 



d«r ldj , (17) 

where r^ is the distance between the unit cells, the term in curly parentheses can now be simplified using the dipole 
approximation to become, 

E E v ij (-i) i+ \-iy+' = e ^^, (is) 

i=l,2j=l,2 e i £ j 



and thus, 



Jmn = E E *w^p-*(^)- (is) 



Hi is the transition dipole moment of a dimer, defined by 



Ml = (20) 



and a is the bond length. Since /ii is independent of the length of the polymer, Eq. (fl9|) reveals how the exciton 
transfer integral varies with the polymer length, as will be described in the next section. However, an alternative 
representation of Eq. (|19[) will now be derived which makes a more direct comparison to the point-dipole approximation 
and introduces the line-dipole approximation. 

Using the expression for the transition dipole moment of the whole chain (from Eq. ([5]), Eq. (fTU|) , and Eq. (HU)) 7 , 
namely 



UN 



K^ViRt), (21) 



where the sum is over unit cells, Eq. ()19|) can be re-expressed as 



2 



4 



This expression for the transfer integral between conjugated polymers justifies the line-dipole approximation as may 
be seen as follows. Using Eq. (|2f p . the transition dipole moment for a polymer of N monomer units may be written 
as a sum of the transition dipole moments for each monomer, namely 

m = J2KRi), (23) 

e 

where 

rtRt)=i*x*(IU)=l*N *V*fly (24) 

Assuming that each monomer of polymer m interacts with each monomer of polymer n via the point-dipole approxi- 
mation, the total exciton integral between the polymers is^, 



4 n =^^ ii{R ti )^n{R tj ). (25) 



Eq. ([221) is thus derived by inserting Eq. $Q into Eq. ((231) . 

Before discussing the evaluation of Eq. (|19|) and its validity, we describe the standing- wave and plane- wave approx- 
imations for the center-of-mass wavefunction, ty(Ri). 



A. Standing waves 



For linear polymers with N equivalent monomers and open boundary conditions the exciton center-of-mass wave- 
function satisfies, 



*(^) = y^ I sin(W, (26) 

where (3 is the pseudo-momentum of the center-of-mass, defined by [3 = irj /(N + 1) with 1 < j < N. This form of 
^(Re) in Eq. (JH]) gives 



m= Mi\/^yCot(/3/2) (27) 



'N/j, 

when N/j 3> 1. We therefore see that in the point-dipole approximation J ~ L. 



B. Plane waves 

In the next section it will be expedient to employ the plane wave approximation for ^(Ri) valid for translationally 
invariant systems, 

= exp(^) 

where K is the Bloch momentum of the center-of-mass, defined by K = 2irj/Nd with —N/2 < j < N/2. This 
approximation allows analytical expressions to be derived for the exciton transfer integral, although it is less applicable 
than the standing wave approximation in linear conjugated polymers and chromophores. 



C. Validity of the Line-dipole approximation 

The validity of the line-dipole approximation has already been discussed elsewhere (see ref^ and references therein). 
When R = d there is a significant error in the line-dipole approximation result. However, for R > 2d the line-dipole 
approximation result is in good agreement with the exact result for short chains, becoming essentially exact for all 
chain lengths when R > 3d. 
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D. Parallel Chains 



An analytical expression can be derived for the transfer integral in the continuum limit assuming the plane wave 
form for the center-of-mass wavefunctions, Eq. (|28|) . Then, 



a\ f L f L (1^3 cos 2 9) , , 
J-J^l I ~ Is J -d Xl dx 2 , (29) 



JO 



where r is the distance between the segments at x\ on chain 1 and xi on chain 2. For parallel chains with D — Do = 
the computation of Eq. (j29|) gives, 



J= (M_) d * ^ 



(30) 



The dependence of the exciton transfer integral on the polymer length is made more explicit by defining J as, 

;=(M) 9 (i>, (3D 
where <?(£), the dimcnsionless exciton transfer integral, is 

g(L) = I I 1 - ; = I (32) 



L 



(1 + ^ 2 ), 



and L — L/R. 

g(L) is plotted versus L for R = 5d in Fig. [Hand Fig. [3] As L -> (or i << i?) g(i) Z/2 and thus J ~ ii? -3 , 
confirming the point-dipole approximation for the whole chain. Conversely, as L — > oo (or L >> R) g(L) — > L _1 and 
thus J ~ (RL)^ 1 . The maximum value of <? is <7 max ~ 0.300, occurring at 



r-\ 1/2 



1.2720- •• (33) 



These idealized analytical results are consistent with the quantum chemistry calculations on explicit polymer systems 
of other authors 2 ^^. 



It is also instructive to express J as, 



where 



./=(§)/.!). (84) 




f(L) = — | 1 - _ | (35) 

'(1 + L*) 

can be regarded as a correction to the point-dipole approximation for the whole chain. As L — > f(L) — > 1, while as 
L —> oo f(L) -> 2L- 2 . /(X) is plotted in Fig. SI 

The analytical results derived in this section have assumed the plane-wave approximation for the exciton center- 
of-mass wavefunction. A more realistic approximation for linear polymer chains, however, is the standing-wave 
approximation, Eq. (|26[) . Analytical expressions cannot be derived using this form of 5*. However, Eq. (JTSJ) is easily 
evaluated numerically. Fig. [5] and Fig. [3] show g(L) in this approximation. For L < 20 the plane and standing wave 
approximations agree rather well. However, the asymptotic behavior is different, decreasing more rapidly with chain 
length for the standing wave approximation, namely g{L) ~ L~ l s , as indicated in Fig. [3j 
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E. Collinear Chains 



For collinear chains with R — 0, Da = L, and D arbitrary, the evaluation of Eq. (|29|) for the plane-wave approxi- 
mation gives, 

J= ~{~D^d) (L + D)(2L + D)' (36) 

or 

where 

g(L) = -x ^ (38) 

(L + 1)(2L + 1) 

and L — L/D. g(L) is plotted versus L/D for D = 5d in Fig. [5] and Fig. [51 As L — + g(L) —> L, while as L — >■ oo 
<?(-£) — > (2L)^ 1 . The maximum value of g is <jf ma x ~ 0.172, occurring at Z max = 1/V2. 

Fig. [5] and Fig. [5] also show the numerical evaluation of the transfer integral using the standing- wave approximation 
for the exciton center-of-mass wavefunction. In this case, as L — + oo g(L) ~ L~ 2 . 

For collinear chains J can also be expressed as a correction to the point-dipole approximation for the whole chain 
as, 



\ D 3 



f(L), (39) 



where 



This function is plotted in Fig. [J] 



f(l) = . (40) 

(L + 1)(2L+1) 



III. DISCUSSION AND CONCLUDING REMARKS 



In this paper the line-dipole approximation for the evaluation of the exciton transfer integral, J, between conjugated 
polymer chains has been rigorously justified for a chain of dimers. Quite generally, it is physically reasonable to apply 
this approximation to a conjugated polymer or linear chromophore with general monomeric units. 

Then, using the line-dipole approximation, as well as the plane-wave approximation for the exciton center-of-mass 
wavefunction, it has been shown analytically that J ~ L when the chain lengths are smaller than the separation 
between them, where L is the polymer length. However, when the chain lengths are larger than their separation, 
J is a decreasing function of chain length, scaling as I -1 for the plane- wave approximation. Scaling relations have 
also been obtained numerically for the more realistic standing-wave approximation for the exciton center-of-mass 
wavefunction, where it is found that for chain lengths larger than their separation J ~ L~ l s or J ~ L~ 2 , for parallel 
or collinear chains, respectively. 

The scaling of J ~ L for short chain lengths is a reflection of the fact that at these lengths scales the point-dipole 
approximation may be applied to the entire chain, implying that J oc (GS|/i|EX) 2 ~ L. Similarly, the scaling of 
J ~ L^ 1 for collinear chains in the plane-wave approximation is easy to understand for chain lengths large compared 
to their separation. In this case the exciton dipoles are uniformally distributed along both chains of length L. The 
double line integral of r~ 3 trivially gives the L~ l scaling. 

The scaling of J with L for long parallel chains is rather counter-intuitive, however, as it implies (as discussed 
below) that the probability of exciton transfer between neighboring parallel chains is a decreasing function of length. 
The physical origin of this behavior can be understood in two ways. First, consider the exciton transfer integral 
for parallel chains when their separation is of the order of the monomer length (d). In this limit J is a periodic 
function of the relative shift variable D^. Fig. [7] shows J for D — (the in-phase configuration), D = a (the out of 
phase configuration), and D — a/2 (the phase-cancelation configuration), illustrating this periodic variation. As the 
separation between parallel chains increases this periodic variation vanishes because of interference effects from longer 
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distance transition densities. It is these interference effects that mean ultimately the exciton transfer integral must 
vanish in the asymptotic limit for chain separations greater than a few monomer lengths. 

A second way to motivate this scaling behavior is to consider the sign of J as two finite-length dipoles are moved 
longitudinally relative to another. As they slide passed each other the sign of J will change once the magic angle 
(6 = cos _1 (l/v / 3)) is reached. Evidently in the limit of infinitely long chains this thought-experiment is nonsensical, 
implying that J = 0. 

The results for the exciton transfer integral have been derived using the assumption that the exciton center-of-mass 
wavefunction, is described either by a standing wave or a plane wave. The standing wave approximation is more 
relevant than the plane wave approximation in linear conjugated polymers or chromophores, although the plane wave 
approximation enables analytical results to be derived. In reality however, disorder or self-trapping will modify the 
exciton center-of-mass wavefunction from both of these idealized functional forms, and under these circumstances it 
is necessary to evaluate Eq. (fH?)) with the explicit form of Such a calculation is beyond the scope of the present 
paper. 

The results presented in this paper have important implications for the photo-physics of conjugated polymers and 
self- assembled molecular systems. First, as already reportedi^^, the Davydov splitting in aggregates is a decreasing 
function of chain length for chain lengths larger than their separation. Second, the Forster transfer rate, k, is reduced 
by a factor of f 2 (L), where f(L) is the correction to the point-dipole approximation for the entire chain, shown in 
Fig. 2] This conclusion may be understood by noting that in the point-dipole approximation, 

1 f°° I D (uj)a>A(oj) , 
k(x ™ -3 duj > ( 41 ) 



R 6 . 

where Iu(ui), the donor emission spectrum, and a^cj), the acceptor absorption spectrum, are both proportional 
to the square of the transition dipole moments. In the line-dipole approximation, therefore, there is a correction 
of f(L) for each of these squared transition dipole moments. This latter result has obvious important deleterious 
consequences for the performance of polymer photovoltaic devices where exciton migration occurs via Forster transfer 
between neighboring conjugated polymers or chromophores. 
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Figure Captions 

FIG 1: The geometry of two conjugated polymers. R is the transverse separation and Do + D is the longitudinal 
shift. Also shown are the sign of the transition densities for the lowest excited exciton, the site labels '1' and '2' in a 
unit cell, and the definitions of a, d, L, and 9. For parallel chains Do — 0, while for collinear chains R — and Dq = L. 

FIG 2: The dimensionless exciton transfer integral, g(R,L) = J R 2 d/2n\ : evaluated for parallel chains (Do = D = 0) 
using the line-dipole approximation with the plane-wave approximation Eq. (|28|) (solid curve) and the standing 
wave approximation Eq. (|26p (dashed curve) for the exciton center-of-mass wavefunction. The inter-chain separation 
R = 5d, where d is the monomer length. 

FIG 3: As for Fig. El 

FIG 4: The correction function, f(L), for the point-dipole approximation. Solid curve: parallel chains where 
L = L/R; dashed curve: collinear chains where L = L/D. 

FIG 5: The dimensionless exciton transfer integral, g(D,L) = J D 2 d/2n\ : evaluated for collinear chains (R = 
and Do = L) using the line-dipole approximation with the plane-wave approximation Eq. (|28p (solid curve) and 
the standing wave approximation Eq. (|26[) (dashed curve) for the exciton center-of-mass wavefunction. The chain 
separation D = 5d, where d is the monomer length. 

FIG 6: As for Fig. [5] 

FIG 7: The dimensionless exciton transfer integral, g(R, L) = JR 2 d/2/ J i 2 , evaluated for parallel chains (Do — 0) using 
transition densities (solid lines) and the line-dipole approximation (dashed lines). The inter-chain separation, R = d. 
The longitudinal shift, D, is D = (circles), D = a/2 (no symbols), and D — a (diamonds). 
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